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The size of higher twist corrections to the spin proton and neutron g\ structure 
functions and their role in determining the polarized parton densities in the nucleon 
is discussed. 



1. Introduction 

The study of the quark-hadron duality phenomena 1 using the present much 
more precise data for the polarized and unpolarized structure functions is 
in a progress 2 . The better understanding of the non-perturbative higher 
twist effects is important for this analysis, especially in the polarized case, 
where the investigations are in the very beginning. 

In this talk we will discuss both, the size of the higher twist effects in 
polarized DIS as well as their role in the determination of the polarized 
parton densities (PPDs) in the nucleon using different approaches of QCD 
fits to the data. 
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2. QCD Treatment of the g± Structure Function 

In QCD the spin structure function gi can be written in the following form 
(Q 2 » A 2 ): 

9i{x,Q 2 ) = gi(x,Q 2 ) LT + g 1 (x,Q 2 ) HT , (1) 

where "LT" denotes the leading twist (r = 2) contribution to gi, while 
"HT" denotes the contribution to g\ arising from QCD operators of higher 
twist, namely r > 3. In Eq. (1) (the nucleon target label N is dropped) 

9i(x, Q 2 )lt = gi(x, Q 2 ) pQ cD + h TMC (x, Q 2 )/Q 2 + 0(M 4 /Q 4 ) , (2) 

where h TMC (x, Q 2 ) are the calculable 3 kinematic target mass corrections, 
which effectively belong to the LT term. gi(x, Q 2 ) p qcd is the well known 
(logarithmic in Q 2 ) pQCD expression and in NLO has the form 

Nf 

.9i(z,Q 2 )pQCD = ^e 2 [(A g +A g -)0(l + ^M^) + ^^AG®^], 

q 

(3) 

where Aq(x, Q 2 ),Aq(x, Q 2 ) and AG(x, Q 2 ) are quark, anti-quark and gluon 
polarized densities in the proton, which evolve in Q 2 according to the spin- 
dependent NLO DGLAP equations. SC(x) q ,G are the NLO spin-dependent 
Wilson coefficient functions and the symbol <S> denotes the usual convolution 
in Bjorken x space. Nf is the number of active flavors. 
InEq. (1) 

9i(x, Q 2 )ht = h(x, Q 2 )/Q 2 + 0(l/Q 4 ) , (4) 

where h(x, Q 2 ) are the dynamical higher twist (t = 3 and r = 4) corrections 
to g\, which are related to multi-parton correlations in the nucleon. The 
latter are non-perturbative effects and cannot be calculated without using 
models. That is why a model independent extraction of the dynamical 
higher twists h(x, Q 2 ) from the experimental data is important not only for 
a better determination of the polarized parton densities but also because 
it would lead to interesting tests of the non-perturbative QCD regime and, 
in particular, of the quark-hadron duality. 

One of the features of polarized DIS is that a lot of the present data 
are in the preasymptotic region (Q 2 ~ 1 — 5 GeV 2 , 4 GeV 2 < W 2 < 
10 GeV 2 ). While in the unpolarized case we can cut the low Q 2 and W 2 
data in order to minimize the less known higher twist effects, it is impossible 
to perform such a procedure for the present data on the spin-dependent 
structure functions without losing too much information. This is especially 
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the case for the HERMES, SLAC and Jefferson Lab experiments. So, to 
confront correctly the QCD predictions with the experimental data and 
to determine the polarized parton densities special attention must be paid 
to the non-perturbative higher twist (powers in l/Q 2 ) corrections to the 
nucleon structure functions. 

3. Higher Twists and Their Role in Determining PPDs 

We have used two approaches to extract the polarized parton densities 
from the world polarized DIS data. According to the first 4 the leading 
twist LO/NLO QCD expressions for the structure functions g\ and F\ have 
been used in order to confront the data on spin asymmetry A\{k, gi/F\) 
and gi/F\. We have shown 5 ' 6 that in this case the extracted from the world 
data 'effective' HT corrections h 9l / Fl (x) to the ratio g\jF\ 

' 9i(x,Q 2 ) 
Fi(x,Q 2 ) 

are negligible and consistent with zero within the errors, i.e. h 9l / Fl (x) « 0, 
when for (#i)lt and (F\)lt their NLO QCD approximations are used. 
(Note that in QCD the unpolarized structure function F\ takes the same 
form as gi in (1), namely F\ = (F\)lt + {F\)ht-) What follows from 
this result is that the higher twist corrections to g\ and F\ approximately 
compensate each other in the ratio g\ jF\ and the NLO PPDs extracted this 
way are less sensitive to higher twist effects. This is not true in the LO 
case (see our discussion in Rcf. 7). The set of polarized parton densities 
extracted this way is referred to as PB(gf LO /F^ LO ). 

According to the second approach 7 , the g\/F\ and A\ data have been 
fitted using phcnomcnological parametrizations of the experimental data 
for the unpolarized structure function F 2 (x, Q 2 ) and the ratio R(x,Q 2 ) of 
the longitudinal to transverse 7^V cross-sections (i.e. F\ is replaced by its 
expression in terms of usually extracted from unpolarized DIS experiments 
F 2 and R). Note that such a procedure is equivalent to a fit to (gi) e xp, but it 
is more consistent than the fit to the gi data themselves actually presented 
by the experimental groups because here the gi data are extracted in the 
same way for all of the data sets. In this case the HT corrections to gi 
cannot be compensated because the HT corrections to Fi(F 2 and R) are 
absorbed in the phenomenological parametrizations of the data on F 2 and 
R. Therefore, to extract correctly the polarized parton densities from the 
gi data, the HT corrections (4) to gi have to be taken into account. So, 



_ 9i(x,Q 2 )lt , h^/^(x) 

<=> TT-, T^TTs H 



F\(x, Q 2 )lt 



Q 2 



(5) 
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according to this approach we have used the following expression for the 
ratio <7i /F\ : 



F 1 N (x,Q*)_ 



g^(x 1 Q 2 ) LT + h N (x)/Q 



F»{x,Q* 



(6) 



exp 



where gi(x, Q 2 )lt (N=p, n, d) is given by the leading twist expression 
(3) in LO/NLO approximation including the target mass corrections. In 
(6) h N (x) are the dynamical r = 3 and r = 4 HT corrections which are 
extracted in a model independent way. In our analysis their Q 2 dependence 
is neglected. It is small and the accuracy of the present data does not allow 
to determine it. The set of PPDs extracted according to this approach is 
referred to as PD(^ T + HT). The details of our recent analysis using the 
present available data on polarized DIS are given in 8 . 

The extracted higher twist corrections to the proton and neutron spin 
structure functions, h p (x) and h n (x), are shown in Fig. 1. As seen from Fig. 
1 the size of the HT corrections is not negligible and their shape depends 
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Figure 1. Higher twist corrections to the proton and neutron g\ structure functions 
extracted from the data on g\ in NLO(MS) QCD approximation for gi(x,Q 2 ) LT - The 
paramctrization (7) of the higher twist values is also shown. 



on the target. In Fig. 1 our previous results on the higher twist corrections 
to gi (before the JLab Hall A data were available) are also presented. As 
seen from Fig. 1, thanks to the very precise JLab Hall A data 9 at large x 
the higher twist corrections to the neutron spin structure function are now 
much better determined in this region. In Fig. 1 our parametrizations of 
the values of higher twists for the proton and neutron targets 
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h p (x) = 0.0465 



h n {x) = 0.0119 



0.1913 




exp[-2((x - 0.2087)/0.2122) 2 



exp[-2{(x - 0.0783)/0.1186) 2 



(7) 



are also shown. These should be helpful in a calculation of the nucleon 
structure function g\ for any x and moderate Q 2 in the experimental region, 
where the higher twist corrections are not negligible. The impact of the 
very recent COMPASS data 13 on the values of higher twist corrections is 
negligible. The new values are in a good agreement with the old ones 8 . 

The values of the higher twist corrections to the proton and neutron g\ 
structure functions extracted in a model independent way from polarized 
DIS data are in agreement with the QCD sum rule estimates 10 as well as 
with the instanton model predictions 11 but disagree with the renormalon 
calculations 12 . About the size of the HT corrections extracted from the 
resonance region see the discussion in the Fantoni's talk at this Workshop. 
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Figure 2. NLO(MS) polarized parton densities PD((;^ T + HT) (solid curves) together 
with their error bands compared to PD(pf LCypNLO) (dashed curves) at Q 2 = 4 GeV 2 . 



In Fig. 2 we compare the NLO(MS) polarized parton densities PD(g^ T + 
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HT) with PD(^ LO /F 1 NLO ). As seen from Fig. 2 the two sets of PPDs are 
very close to each other, especially for u and d quarks. This is a good 
illustration of the fact that a fit to the g\ data taking into account the 
higher twist corrections to g\ (Xdf.nlo = 0-872) is equivalent to a fit of 
the data on ^4i(~ g\/F\ ) and g\/F\ using for the g\ and F\ structure 
functions their NLO leading twist expressions (xdf nlo = 0-874). In other 
words, this fact confirms once more that the HT corrections to g\ and Fi 
approximately cancel in the ratio gi/F\. Nevertheless, we consider that the 
set of the polarized parton densities PD(g]" T + HT) is preferable because 
using them and simultaneously extracted higher twist corrections to g\ , the 
spin structure function gi can be correctly calculated in the preasymptotic 
(Q 2 , W 2 ) region too. 

In conclusion, the higher twist effects in polarized DIS have been stud- 
ied. It was shown that the size of the HT corrections to the spin structure 
function g\ is not negligible and their shape depends on the target. It was 
also demonstrated that their role is important for the correct determination 
of the polarized parton densities. 
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